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ABSTRACT 

The exclusive decay modes B A77 and Bg viv')'!'! shown to have signihcant 
branching ratios of approximately 0.5 x 10“^. This hrst calculation of these modes employs 
a model based on a cascade transition R —>■ I/7 —P77 for estimating the long-distance 
contribution and the process b —> 577 for the short distance one. 
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The investigation of flavor-changing weak radiative transitions of 6-quark has emerged 
during the last decade as a most fruitful held of research, both experimentally and the¬ 
oretically. There are already two measurements by CLEO of transitions in which the 
electromagnetic penguin 6 —sy plays a dominant role, of the inclusive decay transition 
BR{B — >■ Xg'j) = (2.32 ± 0.57 ± 0.35) x 10“*^ |]^ and of an exclusive transition BR{B —> 
77*7) = (4.5±1.5±0.9) X 1O“^0. These measurements conhrm within the existing accuracy 
the Standard Model (SM) prediction, including QCD corrections [Q], for these modes. In par¬ 
ticular, with the inclusion of next to leading order QCD corrections!^ one arrives at the SM 
theoretical prediction BR{B —>• = (3.28 ± 0.33) x 10““^. Consequently, the study of 

these decays is also a good testing ground for new physics^. Concerning the long-distance 
contributions, we remark that various theoretical treatments 0 agree that such contributions 
are small compared to the short-distance electromagnetic penguin b —> sy, most probably 
amounting to less than (5 — 8)% of the rate. 

The two-photon decays of 5-mesons, related to the quark transition b —>■ syy, are also 
of considerable interest]^, |], p. Considering only the contributions of the quark transition 
b —*■ syy, without QCD-corrections, these authors found a fairly sizable rate within SM for 
the two-photon decay of Bg, BR{Bs —> yy) — (1.5 — 3) x 10“^, which should be measurable 
at future 5-machines. 

In a recent paper [jl^, a detailed study of the inclusive rare process 5 —> X^yy has 
been undertaken in terms[|^ of the quark level transition b —>■ syy, both in SM and in two 
versions of the Two Higgs Doublet model. They hnd in SM a branching ratio of ~ 1 x 10“^ 
for decay into hard photons, of energies above lOOMeV each. 

In the present paper we propose a model for calculating the exclusive decays of 5- 
mesons into one pseudoscalar meson (5) and two photons, concentrating on the modes 
Bg —>■ yyy, Bg —> y^y, and B~^’^ —>■ X+’^yy. The calculation includes both short-distance 
and long-distance contributions. Our results reveal two interesting features of these modes: 
hrstly, they appear to contribute a large proportion of the inclusive rate, as estimated in Ref. 
Iim , and secondly, the short- and long-distance contributions are of comparable magnitude. 
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Decays of type B{Bs) —D 77 , involve different features and will be discussed in a separate 
publication]^ . 


We start with the consideration of the b —> 577 transition, which is at the core of our 
treatment. The basic amplitude at the quark level is given Bia by the expression: 


A{b —>• S77) = —i 


OieCp 

71“ 


71 


e^i{ki)eu{k2)u{ps)T^''u{pb), 


( 1 ) 


where 


^ XiTr = Xu{Tr-Tn + MTr-Tn, ( 2 ) 

i=u,c,t 

and Xi = VibV*g. These tensors are divided into the one-particle-irreducible (IPI) and the 
one-particle-reducible (IPR) parts. The IPI part is induced by diagrams with both photons 
coming from the quark-W boson loop, while the IPR diagrams have one photon emitted 
from the external 6 or s quark line. 

For concreteness, we present the details of our approach by treating the Bg —> 777 
decay. Afterwards, we shall refer to three additional modes, Bg —> 7^77 and 
The IPI contribution constituting the short-distance part of the exclusive decay, can be 
calculated directly. The long-distance part of it may be represented by Bg —>■ pcV 77h 
and we estimate this process to be less important. On the other hand, the IPR part needs 
to be performed at the hadronic level and will provide the long-distance part of the exclusive 
transition. For this part, we construct a model which uses vector-meson-dominance (VMD), 
so that the hnal 777 state is realized via the cascade decay Bg —> 07, 0 —> 77. Both these 
intermediate transitions have sizable strength, which is essentially under control (0 —> 77 
is measured and Bg —> 07 is related to the observed B K*'j decay). Accordingly, we 
consider our model to be a reliable tool for the estimation of B{Bg) —>■ P'j'j decays. We also 
made rough estimates for other possible cascade mechanisms, like Bg —>■ B*^ 777 , and 

we hnd this contribution to be signihcantly smaller. 

The IPI part of the transition b 577 with on-shell photons with momenta /ci,fc 2 is 
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given byi, 0 




- k2)i + - 6"«"“A:2^)7„7i 

/Cl • /C 2 


where 72 , = (1 — 75 )/ 2 , = 2 fci ■ k 2 /m^, andjlS 

2 /■! du 


z nil 

< 53 ( 2 :^) = IH-/ —log [1 - 2 :iM(l - m)] . 

Zi Jo u 


( 3 ) 


(4) 


This function ^ 3 (^ 4 ) has been explored in details in [^. 

To calculate the decay amplitude of the IPI part, we need the hadronic matrix element 


{v\siajLb\B,) = ^f+{q^){PB, +Pn)a + ^f-{q^){PBs - Pv)a, 


(5) 


= {Pbs ~ P'nY- For the formfactors /± at zero recoil we use the chiral perturbation 
theory!!^ which gives 

fSs PB ^ Bsv '^ Bs 


f+{{mB, -m^) ) = -f-{{mB, - m^) ) = 


L A+ m„ Vs’ 


( 6 ) 


where A = tub* — mB^- Away from the zero recoil point a monopole behavior will be used: 


2 n r ff . 2 ^^- i .' rnBs -' m ^) / m %, 


f±{q) = f±{{mB^-m^))- y 2 

i — g / 

The result for the short-distance part of 777 is thus given by 

oi^Gp 8 , 


A^^\Bs 777 ) = 


-[A„((53(2:„) - 52,{Zc)) + - 5■i{Zc))\f+{q )e^,uap 


(7) 


( 8 ) 


y27r 9 

•[€^ 62(^1 - (62 ■ kie^ki''k2°‘p^ - ei ■ k2ei^ki’'k2°‘p^)/{ki ■ k2)]. 


We turn now to the reducible part of the amplitude, which we assume to proceed via 
-Bs 07 h 77 ' The amplitude for Bg —>■ 07 is |0 

Gpemb 


A{Bs ^ (t)^{ki,ei)) = 


2 ^ 2 - 


TT^ 


'-VaV:.Cp[Tmie“^^^UvPB.,h. 


(9) 


+r2(0)((e* ■ ti)(PB. ■ ii) - (f* ■ kOipB, ■ ei))l, 

where = 0.65, and Ti(0) = T 2 ( 0 ) = 0.115. This decay is driven by the 6 —»• 57 transition. 


dehned by 


2G 

A{b 57 ) = — -^XtGj^^{mb)Oj, with O 7 = 
V 2 


16ir2 2 


( 10 ) 
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We note that from Eqs. 
comparable to that of B 
The amplitude for 0 


(9) and (10) one predicts]^ a branching ratio for Bg 
> K*'y, which is indeed expected. 

->• 777 is 


07 


A{(j) 77 (/c2 , € 2 )) = (11) 

where |c 0 ,,.y| = 0.21GeV“^ is determined from the partial decay width. Thus the IPR part 
of the amplitude turns out to be in the VMD model 

A^^^iBg ^ 777 ) = + T 2 ( 0 )(ei„PB. ■ k^ - k.^Ps^ ■ e^] 

(^2-yPr,sk2p + im^V^ 

+ (ki > ^ 25^1 ^ 2 )- ( 12 ) 


This treatment of the IPR amplitude is consistent with the decomposition theorem^, |T^ . 

We turn now to extend our calculation to include also the Bg —>■ 7^77 mode. For this 
purpose, if we use the nonet symmetry between the octet and the singlet light pseudoscalars, 
we need to consider the p — p' mixing, which is dehned by 


p = pscosOp — 7osin6*p (13) 

p' = pgsinOp + pocos9p, 


and the factor \j2/3 in Eq. ( 6 ) should be repaced by ^J2J3cos9p + ^Jl/ Ssin^p for the Bg ^ p 
transition, and by y^2/ 3sin0p — y^l / 3cos9p for the Bg p' transition. Numerically 9p ~ —20° 
will be used. Now, the physical masses of p or p' will be used in Eq. ( 6 ). In this way, Eq. 
( 8 ) is extended to cover the decays to p and 7 ', with the appropriate replacements in ( 6 ) and 

( 7 ). 


For the long-distance part of Bg p'pp, we use again the nonet symmetry, under which 
the coupling for 0 —>• p'p is related to by||T7[ 


\c<i,p'-f/c^p^y\ = {cos9p - V2sm9p)/{sin9p V2cos9p), 


(14) 
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which gives \ctf,n'^\ = 0.30GeV“^ if 9p = —20° is used. A completely similar calculation is 
carried out for B K'J'J, in which case the long-distance part is due to —>■ + 

Here, the required branching ratios are known for both transitions [|T^] . 

In presenting our results, we recognize that unfortunately we cannot £x the relative 
phase between the IPR and the IPI amplitudes. This makes the phase of the interference 
effects in the decay rates undetermined. The reason can be seen from the amplitudes in 
(9) and in (12). There are too many Lorentz structures involved, which make the analogue 
Argand plot analysis]^ invalid. Thus we will present the results for the IPI, the IPR and the 
interference contributions separately . The total decay rates are then Tipj + Tipji±Tinter- As 
a result, we also refrain from presenting differential decay distributions in various variables. 

We follow now the authors of Ref. |l^ by imposing several cut conditions on the 
decays considered, which were enforced on the quark level transition b —>■ s'j'j in Ref. [jT0|. 
First, hard photons with energy larger than lOOMeV are chosen. Second, the invariant mass 
squared for any two hnal particles is demanded to be larger than cm^ or with c =0.01 

or 0.02. Third, all the angles between two of the hnal particles are taken to be larger than 
20°. In addition, we wish to exclude from the decay the region of and iPy which is 

close to the resonance peak. Thus, we demand that the invariant mass of any r]'y,r]''y pair 
deviates from the 0 resonance by more than ±50(or 100)MeV, while for K'y the requirement 
extends to 300(or 500)MeV (the width of 0 is 4.43MeV and of K* is 50MeV). We also 
used a modihed form of the new requirement by demanding that the invariant mass of any 
v{v')l be larger than l.lGeV or 1.2GeV, and of K'y be larger than 1.2GeV or 1.4GeV. 
Using the same mass parameters and GKM matrix elements as in [^, and = 0.13GeV, 
fri' = O.llGeV, fx = O.lOGeVpPI, gB*Bsv — 0-^) fhe results are given in Table 1. In the 
calculation we approximate the strengths of various vertices, {BK*'y), {(pg'j), etc. with their 
values on the mass-shell. 

We have studied the sensitivity of our results on the cut conditions. To do this, we £x 
the cut condition to assure the offshellness of the 0 resonance and then modify the other 
conditions. In Table I we present numerical results for representative cuts. First, weak 
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dependence is found for the hard-photon requirement and c; with or without these two 
conditions, the decay rates vary within 2 % for both Bg 1777 and Bg —> 7 ' 77 . Second, the 
dependence on the angular cut is moderate; even using a cut as large as 20 °, the reductions 
for the IPI and the IPR rates are within 5%, while no visible change in the interference effect 
exists. Indeed, the interference effect is almost stable against all different cut conditions. 

The insensitivity on the cut conditions can be understood from two aspects. First, 


unlike the quark level transition used in |T^ where the IPR amplitude is quite singular, 
these singularities are absent in the present model. Second, the new requirement of the 
offshellness of the intermediate vector meson eliminates quite a portion of the phase space 
for small angle between the P and 7 . Although it is convenient to use also the cut conditions 


of in the experimental analyses of the signals, these conditions are less important on the 
theoretical side in our case. 

The only sizable dependence of the decay rates is on the cut on the vector meson 
resonance. The strongest dependence is for the IPR rates, which vary almost 20% for 
Bg —>■ 777 under the different conditions enforced. For Bg v'll fhis dependence is 
much weaker, since the strong interference of the two possible 0 conhgurations occurs away 
from and dominates over the resonant region. Numerically, the typical decay rates are 
(1.58 ± 0.18) X 10“^°GeV for Bg —>• 777 and (2.42 ± 0.48) x 10“^°GeV for Bg —>• ^' 77 , if 
~ — lOOMeV is required. The errors given here, as well as those given for 

B —>■ 7 ^ 77 , are due to the uncertainties in the phases of the interference terms. Taking 
T{Bg) = 1.61 X 10“^^s, the corresponding branching ratios are (0.39 ± 0.04) x 10“^ for 
Bg —>• 777 and (0.59 ± 0 . 12 ) x 10 “^ for Bg —>• ^' 77 , respectively. Also, the averaged open 
angles between the two photons for these two channels are all around 140°, which is stable 
against the phases of the interference effects and is comparable to the value ~ 135° in [lU . 
For B —^ A" 77 , again, the only sizable dependence of the decay rates is on the cut on the 
vector meson resonance. The decay rates are (2.87 ± 0.39) x 10“^°Gev for B^ —> 
and (2.04 ±0.26) x 10“^°Gev for B^ —^ if | — "^k* \ > 300MeV is taken. They 

correspond to branching ratios of (0.68±0.09) x 10“^ for B^ —^ and (0.50±0.06) x 10“^ 
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for The averaged open angle between the two 7 ’s is also around 140°. 

To summarize, our results reveal interesting features of the exclusive decays Bg —>■ 
vivh'l- Firstly, inspection of Table I indicates that the short-distance(lPI) and the long- 
distance(lPR) components play a comparable role in these decays, which is a new feature 
in the domain of B radiative decays. 

Secondly, the calculated exclusive modes appear to constitute a very sizable portion of 
the respective inclusive decays. For example, in [|T^ it was found BR{B —>■ Xg'j'j) ~ 1 x 10“^, 
to which our results for B —>■ K'j'y are compared. We remark, however, that unlike the 
B —> Xg'y transition where the dominant mechanism is the quark level b ^ s'y transition, in 
the case of i? —> Xg'y'y the diagram with one photon attached to the spectator quark also 
contributes. This mechanism is accounted by our treatment but has not been included in 
0 , which might affect the ratios of the calculated exclusive and inclusive decays. 

We would like to acknowledge discussions with Prof. G. Eilam. This research is sup¬ 
ported in part by Grant 5421-3-96 from the Ministry of Science and the Arts of Israel. The 
research of P.S. has also been supported in part by the Fund for Promotion of Research at 
the Technion. 


References 

[ 1 ] M.S. Alam et ah (GLEO Gollaboration), Phys.Rev.Lett.74, 2885 (1995). 

[2] R. Ammar et ah (GLEO Gollaboration), Phys. Rev. Lett. 71, 674 (1993). 

[3] S. Bertolini, F. Borzumati and A. Masiero, Phys. Rev. Lett. 59, 180 (1987); N.G. 
Deshpande, P. Lo, J. Trampetic, G. Eilam and P. Singer, Phys. Rev. Lett. 59, 183 
(1987); R. Grigjanis et ah, Phys. Lett. B213, 355 (1988); B286, 413 (E) (1992); B. 
Grinstein, R. Springer and M. B. Wise, Phys. Lett. 202B, 138 (1988); G. Gella et ah, 
Phys. Lett. B325, 227 (1994). 





[4] M. Ciuchini et al., Nucl. Phys. B415, 403 (1994); M. Misiak and M. Miinz, Phys.Lett. 
B344, 308 (1995). A.J. Buras et al., Nucl. Phys. B424, 374 (1994); K. Chetyrkin, M. 


Misiak and M. Miinz, [hep-ph/9612313 


[5] J.L. Hewett, Proc. of the “1994 International Workshop on B Physics”, Nagoya, Japan, 
Eds. A.I. Sanda and S. Suzuki, World Scientihc, 1995. 


[6] E. Golowich and S. Pakvasa, Phys. Rev. D51, 1215 (1995); G. Ricciardi, Phys. Lett. 
B355, 313 (1995); N. G. Deshpande, X.-G. He and J. Trampetic, Phys. Lett. B367, 362 
(1996); G. Eilam et al., Phys. Rev. D53, 3629 (1996). 

[7] G.-L. Lin, J. Liu and Y.-P. Yao, Phys. Rev. D42, 2314 (1990). 

[8] H. Simma and D. Wyler, Nucl. Phys. B344, 283 (1990); T.M. Aliev and G. Turan, Phys. 
Rev. D48, 1176 (1993). 

[9] S. Herrlich and J. Kalinowski, Nucl. Phys. B381, 501 (1992). 

[10] L. Reina, G. Ricciardi and A. Soni, Phys. Lett. B396, 231 (1997). 

[11] 1. 1. Bigi, Proc. of the “1994 International Workshop on B Physics”, Nagoya, Japan, P. 
391, Eds. A.I. Sanda and S. Suzuki, World Scientihc, 1995. 


[12] P. Singer and D.-X. Zhang (unpublished). 

[13] M.K. Gaillard and B.W. Lee, Phys. Rev. DIO, 897 (1974); E. Ma and A. Pramunida, 
Phys. Rev. D24, 2476 (1981). 

[14] M. B. Wise, Phys. Rev. D45 (1992) 2188. 

[15] P. Singer, Phys. Rev. D49, 7 (1994). 

[16] F. E. Low, Phys. Rev. 110, 974 (1958). 

[17] P. Ko, J. Lee and H. S. Song, Phys. Lett. B366, 287 (1996) . 


9 



[18] R. M. Barnett et al. (Particle Data Group), Phys. Rev. D54, 1 (1996). 

[19] P.J. O’Donnell, and H.K.K. Tung, Phys. Rev. D43, 2067 (1991). 

[20] H. Aihara et al. (TPC /27 Collaboration), Phys. Rev. Lett. 64, 172 (1990). 


10 



TABLE I. Decay rate for Bs Vlli A °77 and A ^77 (in 10 ^'^GeV). 




no angular cut 



all angles > 20° 


Res. Cut 

IPI 

IPR 

Inter. 

IPI 

IPR 

Inter. 

Bs mi 



- m0 > 50MeV 

0.9S 

0.83 

0.18 

0.92 

0.79 

0.18 

> lOOMeV 

0.92 

0.66 

0.18 

0.90 

0.64 

0.18 

^rn^-y > 1.2GeV 

0.87 

0.S6 

0.17 

0.87 

O.SS 

0.17 

Bs v’m 



-m^\> 50MeV 

1.77 

0.66 

0.48 

1.72 

0.63 

0.46 

-m^\> lOOMeV 

1.77 

0.6S 

0.48 

1.72 

0.62 

0.46 

y/^ > 1.2GeV 

1.77 

0.64 

0.47 

1.72 

0.61 

0.46 

R° ^ A°77 



\\J'ni\^ — rriK* > SOOMeV 

1.40 

1.48 

0.39 

1.40 

1.47 

0.39 

Y^W'X7 “ mK* \ > SOOMeV 

1.31 

1.30 

0.37 

1.31 

1.29 

0.37 

-JrnK^ > 1.4GeV 

1.24 

1.19 

0.37 

1.24 

1.19 

0.37 

^ A±77 



\yj'ni\^ — rriK* > SOOMeV 

1.40 

0.64 

0.26 

1.40 

0.64 

0.26 

Y^W'X7 “ rnK* \ > SOOMeV 

1.31 

0.S6 

0.2S 

1.31 

0.S6 

0.2S 

> 1.4GeV 

1.24 

0.S2 

0.24 

1.24 

O.Sl 

0.24 
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